Fission barriers in neutron-proton isospin plane for heavy neutron-rich nuclei 
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We discuss the sensitivity of fission barrier for heavy neutron-rich nuclei to fission paths in the two 
dimensional neutron-proton quadrupole plane. To this end, we use the constrained Skyrme-Hartree- 
Fock + BCS method, and examine the difference of fission barriers obtained with three constraining 
operators, that is, the neutron, proton, and mass quadrupole operators. We investigate U, ' U, 
and 266 U, that is relevant to r-process nucleosynthesis. We find that the fission barrier heights are 
almost the same among the three constraining operators even for neutron-rich nuclei, indicating 
that the usual way to calculate fission barriers with the mass quadrupole operator is well justified. 
We also discuss the difference between proton and neutron deformation parameters along the fission 
paths. 

PACS numbers: 25.85.-w, 26.30.Hj, 26.30. +k, 21.60.Jz 



I. INTRODUCTION 

Fission plays a decisive role in determining the stabil- 
ity of heavy nuclei, where the Coulomb energy competes 
with the nuclear surface energy. Typical examples are 
superheavy elements (SHE). Although the fission barrier 
disappears in SHE in the liquid drop model, the nuclear 
shell effect leads to a relatively high fission barrier and 
eventually stabilizes SHE. The experimental efforts have 
been continued in many facilities to synthesize SHE using 
heavy- ion fusion reactions 

It has been well recognised that fission plays an impor- 
tant role also in the context of nuclear astrophysics 0,03], 
but systematic investigations on its role in r-process nu- 
cleosynthesis have started only in recent years 0, [H, @|. 
The r-process is one of the most promising candidates 
that synthesize the elements heavier than iron (Fe). In 
this model, nuclei capture a number of neutrons via suc- 
cessive (n, 7) reactions in a highly neutron-rich environ- 
ment e.g., neutron stars. As a consequence, the r-process 
path passes through the neutron-rich side in nuclear chart 
which cannot be reached experimentally at this moment. 
Heavy neutron-rich nuclei produced by the r-process may 
decay by spontaneous fission, neutron-induced fission, or 
beta-delayed fission 0, 0, H, @|. The neutrino-induced 
fission might also play a role if the neutrino flux is signif- 
icant ■ In order to construct a reliable r-process model 
with fission, it is urged to calculate systematically fission 
barriers of many neutron-rich nuclei. 

Theoretically, fission barriers can be calculated using 
either the macroscopic-microscopic model || or micro- 
scopic mean-field models 0, [l(J EH- ^ n the latter ap- 
proach, one selects a few important degrees of freedom for 
fission, such as quadrupole or higher multipole moments, 
and draws a fission energy surface using the constrained 
Hartree-Fock method with the corresponding constrain- 
ing operators. The total energy is minimized with respect 
to all the other degrees of freedom than those considered 
explicitly in the calculation. In this sense, the mean-field 
approach provides an adiabatic potential energy surface 
for the case where the selected degrees of freedom are 



much slower than the other degrees of freedom so that 
they adiabatically follow the motion of the former at ev- 
ery instant. 

Usually, one takes a mass (i.e., proton+neutron) 
quadrupole moment as one of the most important de- 
grees of freedom. This implicitly assumes either that 
the isoscalar motion is much slower than the isovector 
motion or that the isoscalar and isovector motions are 
decoupled. For fission of neutron-rich nuclei, however, it 
is not obvious whether this assumption is justified, and it 
may be more natural that the shape degrees of freedom 
for neutron and proton are treated separately. In fact, a 
two dimensional energy surface spanned by proton and 
neutron deformations has been drawn recently for light 
neutron-rich nuclei, such as 16 C and 20 O [H [ij, El • The 
difference in neutron and proton deformation parameters 
along a fission path for actinide nuclei has also been in- 
vestigated in Refs. flR fl6|. 



In this paper, we examine the sensitivity of fission bar- 
riers for neutron-rich nuclei to a choice of constraining 
operator in the isospin plane. To this end, we use the con- 
strained Skyrme-Hartree-Fock + BCS method, and con- 
sider the mass, proton, and neutron quadrupole moments 
as the constraining operators. The constrained Hartree- 
Fock method with the proton (neutron) constraint pro- 
vides an adiabatic energy surface for the case where the 
proton (neutron) is much slower than neutron (proton). 

The paper is organised as follows. In Sec. II, we use 
a schematic model and illustrate an example in which 
the constrained-Hartree-Fock method with the three dif- 
ferent constraints give significantly different results from 
each other. In Sec. Ill, we briefly summarize the the- 
oretical framework for constrained Skyrme-Hartree-Fock 
method. Sec. IV presents the results for the fission bar- 
rier, single particle levels, and the proton and neutron de- 
formations along the fission paths for uranium isotopes. 
We then summarize the paper in Sec. V. 
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FIG. 1: Schematic pictures of the constrained Hartree-Fock 
method. The left panel is for the total constraint, in which 
the energy is minimized along the lines of Q p + Q n = const. 
shown by the dashed line. The right panel is for the proton 
constraint, in which the energy is minimized along the lines 
of Q p — const. The corresponding paths are shown by the 
thick solid lines in both the figures. 




Qp+Qn 



FIG. 2: The energy along the paths shown in Fig. \T\ The 
solid and the dotted lines are for the total and the proton 
constraints, respectively. 



II. SCHEMATIC MODEL 

In the constrained^ Hartree-Fock method with a con- 
straining operator O, one minimizes the expectation 
value of 

H' = H-\6, (1) 

where H is the Hamiltonian of the system. The Lagrange 
multiplier A is determined so that the expectation value 
of O becomes a given value Oo . The O can be any one- 
body operator, but usually the mass quadrupole opera- 
tor, Q = Q p + Q n , where Q p and Q n are the quadrupole 
operators for proton and neutron, respectively, is consid- 
ered as one of the constraining operators O when one 
studies fission barriers. The aim of this paper is to com- 
pare such fission barriers with those obtained by using 
O = Q p or Q n . We call the latter scheme "proton (or 
neutron) constraint" while the former "total constraint". 

Before we perform self-consistent calculations, we 
would like to illustrate a possible difference among the 
three schemes for the constrained Hartree-Fock using a 
schematic model. Suppose that we have an energy sur- 
face shown in Fig. 1 in the two-dimensional plane of 
proton and neutron quadrupole moments, Q p and Q n . 
In the total constraint scheme, the energy minimum is 
searched along the line (Q) = (Q p ) + (Q n ) = const., 
which is shown by the dashed lines in the left panel. The 
resultant path is denoted by the thick solid line, and the 
energy variation along this path is shown in Fig. 2 by 
the solid line. In the case of proton constraint, on the 
other hand, the energy minimum is searched along the 
dashed lines in the right panel of Fig. 1, which correspond 
to (Q p ) = const. The path and the energy are shown by 
the thick solid line in Fig. 1 and the dotted line in Fig. 
2, respectively. The energy is plotted as a function of the 
total quadrupole moment along the path. Those of the 
neutron constraint are obtained in a similar way. 

We see clearly that the two paths obtained with the 
different constraining operators deviate significantly from 



each other. The energy is also different as a function 
of the total quadrupole moment, although the absolute 
minimum can be obtained irrespective to the choice of 
the schemes. The ambiguity arises when the number of 
degree of freedom is reduced from two to one. 

Notice that the differences among the schemes will be 
small if the energy surface is much steeper along the line 
of (Q p ) + (Q n ) — const. In the next section, we will 
investigate how much the fission barriers are changed for 
realistic nuclei depending on which scheme one employs 
to minimize the energy. 



III. NUMERICAL DETAILS 



In order to calculate fission barriers for realistic nuclei, 
we use the Skyrme-Hartree-Fock+BCS method [TtJ (see 
Ref. P3] for a recent review). In this method, the ex- 
pectation value of the Hamiltonian H is given in terms 
of an energy functional as 



E= drH{r) (2) 
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with 

H(r) = £-t(t) 
2m 



rameters are given by 
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(3) 



Here, p q (r), r q (r), and J q {r) are the nucleon density, the 
kinetic energy density, and the spin density, respectively, 
which are defined as, 

Pl( r ) = V i l^i( r '°''?)| a ' 

J q (r) = (-<) 51 «i #(* , i £r >«)(v&(r > ff / ,gr) x (<r|oV; 



(4) 



In these equations, g denotes the isospin (<?=p or n), ^ 
is the single-particle wave function, and vf is the occu- 
pation probability estimated in the BCS approximation. 
H.c(f) in Eq. ([3]) is the Coulomb energy term, while p, 
t, and J are the total (proton+neutron) densities. 
In this paper, we use the quadrupole operator 
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(5) 



as a constraining operator in Eq. (U). For simplicity, 
we assume the reflection and axially symmetric nuclear 
shapes, although the mass asymmetry sometimes plays 
an important role in describing nuclear fission. From the 
expectation value of the quadrupole operator, we calcu- 
late the total deformation parameter as [19] 
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(Qt), 



(6) 



where i?o is the nuclear radius parameter given by i?o = 
1.1 A 1 / 3 (fm). The proton and neutron deformation pa- 
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(7) 



respectively. 

In the actual numerical calculations shown in the next 
section, we use the computer code SKYAX [13, El- This 
code solves the Skyrme Hartree-Fock equations in the co- 
ordinate space with the reflection and axially symmetries. 
We use the mesh size of 0.6 fm. The pairing correlation 
is taken into account in the BCS approximation. In this 
paper, we use the delta force 



^pair 



{r,r') = -V S(r-r'), 



(«) 



for the pairing interaction. We expect that our conclu- 
sion is qualitatively the same even if we use a density- 
dependent delta interaction. In the code, the smooth 
cut-off function 
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1 + exp ((e Q — A — AE) ///)' 



(9) 



is introduced for the pairing active space. Here, A is the 
Fermi energy, and AE is determined so that 



1.65iV 2/3 , 



(10) 



with p — AE/1Q, N q being the number of particle for 
proton (q — p) or neutron (q = n). We use the strength 
parameter of Vq = 279.082 MeV- fm 3 for proton and 
258.962 MeV ■ fm 3 for neutron pairings [Hj]. 



IV. RESULTS 

We now present the results of constrained Hartree-Fock 
calculation for the fission barriers of 220 < 236 ,266jj nuc i e i 
The 236 U is on the /^-stability line, while 220 U and 266 U 
are proton-rich and neutron-rich nuclei, respectively. No- 
tice that 266 U is relevant to r-process nucleosynthesis [12] ■ 
We adopt the parameter set SLy4 [23| for the Skyrme 
functional. 

Figures [3] and |4] show the fission barriers for the 236 U 
and 266 U, respectively, as a function of the total deforma- 
tion parameter (3t. The upper panels are obtained with 
the total constraint, while the lower panels show the dif- 
ference of the fission barrier obtained with the proton 
constraint from that with the total constraint (the dashed 
line) . A similar quantity for the neutron constraint is also 
shown in the lower panels by the dotted line. The dif- 
ferences are much smaller than the fission barrier height, 
and the fission barriers obtained with the three schemes 
are almost indistinguishable in the scale shown in the 
figure. We have calculated for other even-even uranium 
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FIG. 3: The fission barrier for U as a function of deforma- 
tion parameter (3 t obtained with the total constraint scheme 
(the upper panel). The dashed line in the lower panel shows 
the difference of the fission barrier obtained with the pro- 
ton constraint and that with the total constraint, while the 
dotted line denotes the difference obtained with the neutron 
constraint. 



FIG. 5: Proton single particle levels near the Fermi energy 
for 236 U (the top panel) and 266 U (the bottom panel) as a 
function of the total deformation parameter. The solid and 
the dashed lines are the results for the total and the proton 
constraints, respectively. 
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FIG. 4: Same as Fig. El but for 266 U. 



isotopes from 220 U to 276 U, and confirmed that the three 
schemes lead to almost the same fission barriers for all of 
these nuclei. 

Let us next discuss single-particle levels. Fig. [5] shows 
the proton single-particle energies near the Fermi energy 
as a function of the total deformation parameter. The 
solid and the dashed lines show the results of the total 
and the proton constraints, respectively. We see that the 
single-particle energies are similar to each other between 
the total and proton constraints, although the difference 
is not negligible at large deformations. We have found 
that the tendency is similar also for the neutron con- 
straints, although the deviation is smaller as compared 
to the proton constraint. We have also found that the 
conclusion remains the same also for the neutron single- 
particle energies. 

The difference of deformation parameters for proton 
and neutron along the fission paths is shown in Fig. [6] 
f or 220,236,266tj Although the difference among the three 
curves is now more visible than in the fission barriers, 
the results with the three schemes are similar to each 
other, indicating that the fission path is not sensitive to 
the constraining operator in the isospin space. The (3 P — 
Pn is not a monotonic function of the total deformation 
parameter fit, but on average it increases with /3t for 220 TJ 
and 236 U. Although the average value of (3 P — (3 n appears 
to be zero even for large values of (3 t for the neutron-rich 
nucleus 266 U, this might be an artifact of using the same 
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FIG. 6: Difference between the proton and the neutron defor- 
mation parameters, (3 P — /3 n , for 220 > 236 . 266 u nuclei as a func- 
tion of the total deformation parameter. The solid, dashed, 
and dotted lines are the results for the total, proton, and 
neutron constraints, respectively. 



radius parameter Rq between neutron and proton in Eq. 
([7]) to calculate the deformation parameters (but see the 
discussion below). 

Figure [7| shows the density difference of proton and 
neutron, p p — p n for 236 TJ and 266 1J obtained with the 
total constraint for (3 t = 2.0. It is plotted in the two- 
dimensional plane of (r, z), where the density has the 
axial symmetric shape around the z axis. One can notice 
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FIG. 7: (Color online) Difference of the proton and the neu- 
tron density distributions, p p — p n , for 236 U (the upper panel) 
and 266 U (the lower panel) at (3t = 2.0 obtained with the to- 
tal constraint. The densities are axial symmetric around the 
z-axis. 



that the difference between the proton and the neutron 
densities is larger in 236 U than in the neutron- rich nu- 
cleus 266 U. This is consistent with the difference in the 
deformation parameter shown in Fig. [6J 



V. CONCLUSION 

We have used the constrained Skyrme-Hartree-Fock + 
BCS method with a quadrupole constraint in order to 
calculate the fission barriers of neutron-rich uranium nu- 
clei with astrophysical interests. In particular, we carried 
out the calculations with the proton, neutron, and mass 
(total) quadrupole operators as the constraining opera- 
tors. We have found that the fission barriers are almost 
independent of the constraining operators in the neutron- 
proton isospin space. We have discussed this behaviour 
using a schematic model, and suggested that the poten- 
tial energy surface is steep along the isovector degree of 
freedom. We have also found that the single-particle lev- 
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els as well as the deformation parameters along the fission 
paths do not depend much on the constraining operators. 

Our calculations indicate that the proton and the neu- 
tron deformations differ from each other even for the 
nucleus on the /3-stability line, 236 U, and the difference 
increases as the total deformation parameter becomes 
large. 

In the study of fission barriers based on the mean-field 
approaches, one usually uses mass multipole moments 
as constraining operators. Our results suggest that this 
approach is justified even for neutron-rich nuclei, where 
one might expect that the proton or neutron constraint 
is better. Of course, it is always desirable to treat pro- 
ton and neutron separately for neutron-rich nuclei, and 
draw a two-dimensional fission energy surface. However, 
it is rather demanding to do so if one has to take into 
account explicitly many multipole moments, including si- 
multaneously the quadrupole and octupole moments, or 



even higher multipole moments. Even in that case, our 
study clearly indicates that one can reduce the number 
of degree of freedom by introducing the mass multipole 
moments, rather than treating proton and neutron sepa- 
rately. 
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